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Chapter 1 
Introduction 


Several techniques are available for the determination of constitutive mate- 
rial parameters. The technique of choice is dependent upon sample geome- 
try, frequency range of interest and the relative range of parameter values. 
Some techniques involve the measurement of reflection and transmission 
of RF energy which then can be related to the properties of the geome- 
try samples. Naturally, rigorous relationships have to exist between the 
measured quantities and the material parameters. The parameters which 
describe the materials electrically are the quantities permittivity ( e ) and 
permeability (/x) which relate the electric and magnetic field intensities to 
their flux densities. Mathematically, this relationship is 


D = eE 


( 1 . 1 ) 


and 


B = /xH. (1.2) 

One common apparatus used for the acquisition of the specular reflec- 
tion performance of lossy material samples is the NRL arch. The sample 
is usually a planar conductive panal coated with some material. Attempts 
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for parameter extraction are usually limited to just one parameter for a 
lossy material (usually non- magnetic). Two complex parameters have to 
be extracted if the material is magnetic. 

An accepted measurement technique for constitutive parameter deter- 
mination involves the reflection and transmission measurement of samples 
in an airline or waveguide fixture. This approach, unlike the NRL arch, has 
a known rigorous relationship between the measurable reflected-transmitted 
fields and the constitutive parameters. The practice is to completely fill the 
interior cross section of the fixture with the material for a given length. Sim- 
ple expressions exists to obtain the desired parameters. Potential sources 
of error aside from related hardware limitations are sample dependent. The 
sample uniformity, geometry and position in the fixture can limit the accu- 
racy of the desired parameter values. 

Potential accuracy improvements in the parameter values are possible if 
alternate sample geometries are used. Requiring the sample to completely 
fill the cross section of the fixture requires a good fit between the sample 
and fixture walls. If the fit is not adequate the measured quantities do not 
correspond to the anticipated geometry. Introducing a known gap between 
the sample and a fixture wall may make a less sensitive requirement in 
the actual dimensions for the sample. This gain is offset by the increased 
complexity of the reflected and transmitted field solutions. 

This report discusses and contains results for NRL arch measurements 
and waveguide fixture performances for completely and partially filled waveg- 
uides. 
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Chapter 2 
NRL Arch 


The NRL arch is simply a hardware configuration to acquire the reflection 
performance of material coatings. The term “reflection” refers to scatter- 
ing between the transmit and receive horns which satisfies Snell’s Law of 
Reflection (angle of incidence equals angle of reflection) from the surface of 
the material. This restriction is imposed for the following material to allow 
a description as accurate as possible of the scattered field without a great 
amount of complexity. 

The scattering from a two dimensional geometry as shown in Figure 
2.1 consists of three dominant mechanisms, i.e., two edge diffracted contri- 
butions and one specular contribution when viewed from a high frequency 
viewpoint (the dimension “w” is much larger than the incident wavelength). 
The direct coupling between the transmit and receive horns is ignored here. 
This field decomposition fails when the strip is electrically small. As the 
distance between the transmit-receive horns and the strip become large, 
it is easily shown that the reflected field can be represented by the edge 
diffracted fields alone. Large refers to the distance required to maintain a 
constant phase across the surface of the strip. An accepted distance, R, to 
achieve plane wave illumination (constant phase) is 
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( 2 . 1 ) 


R > 2 


W 2 
X ' 



w 


Figure 2.1: Near field scattering from a strip. 


A detailed study of general bistatic scattering from a perfectly conduct- 
ing strip has been done [1] which analytically and rigorously demonstrates 
the above described behavior. The significance of this analysis is that it 
demonstrates that the scattering from a planar surface which satisfies the 
reflection law can be represented as an edge diffracted field. However, the 
edge diffracted fields for plane wave incidence in the far field does simulate 
a reflected field when reflection law is satisfied. 

The two dimensional strip can be extended into a three dimensional 
plate case when the other orthogonal plate dimension is electrically large. 
The transformation of two to three dimensional fields is given by 


E 3d = 


( 2 . 2 ) 
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and geometrically illustrated in Figure 2.2. 

The principle plane physical optics (PO) approximation for the specular 
scattered planar plate of dimensions w by L with parallel polarization is 


n wLe~’ kr 
E‘ = jR. cos 8— 


(2.3) 


and with perpendicular polarization it is 

n wLe~ ikr 
E' = -jRh cos 0 - . 

A T 


(2.4) 


The polarization designations are with respect to the electric field direction 
and the incident plane of incidence. The quantities R t< h are the Fresnel re- 
flection coefficients for an infinite conducting planar surface with an uniform 
coating of material (thickness being t). The Fresnel reflection coefficients 

are 


jrjr cos 8 tan(fct cos 0*) — cos 6 |t 
jrj r cos 8 tan (kt cos 8 l ) + cos 0‘ 


and 


j rjr cos 8 1 tan(fet cos P) — cos 0 
jrj r cos 8 l tan(fc< cos 8 l ) + cos 8 


(2.5) 


( 2 . 6 ) 


where Tj r = \j^ i ^ = y/J^T^r^o und k„ sin 8 — k sin 8 . 

A more accurate approximation to the specular scattered field is giveii 
by GTD [2,3] based upon edge diffraction through accounting for the first 
order change in the reflection shadow boundary for material coating. The 
scattered field for parallel polarization is 

E * = — — ( — jR f kw cos8) L — — (2-7) 

2n Vcos 8 ) t 
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and for perpendicular polarization it is 


= jRhkw cos ") Lt ir- < 2 - 8 > 

The first term in these expressions is a result of the discontinuity of the 
incident field. 

Conceptually, parameter extraction for both e and fj. is possible when 
measurements for two different incident angles are performed. The two in- 
cident angles would have to be sufficiently different to insure the measure- 
ments to have their own unique character. Attempts have been made based 
upon PO approximations which have yielded poor results [4]. Such results 
would be valid only when the second term is the parenthesis of Equations 
(2.7) and (2.8) is dominant. The GTD approximation would potentially do 
better but the scattered field representation may not be rigorous enough if 
a surface wave is induced in the coating. 

Arch measurements were performed for normal incidence upon planar 
coated metallic plates. The arch was configured with a single 2-18 GHz 
AEL horn 70” from the plates. Better performance can be obtained if two 
horns are used instead of one due to the input mismatch of broadband 
horns. The horn to horn return is up to 20 dB less than the mismatch 
return for a single horn at the higher frequencies. A septum is suggested to 
help minimize the coupling between the horns. A high gain horn will also 
provide greater sensitivity due to lower sidelobe levels to sense the nearby 
environment. 

Figures 2.3 and 2.4 show the measured reflection coefficients for con- 
ducting plate sizes of 5” square and 12” square, respectively, when coated 
with Emerson and Cuming SF-9 material. The average midband material 
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parameters are e r = (7.8,-j.l7) and /x r = (1.42, -j.8) with a thickness of 
.082”. The plots also contain the calculated Fresnel reflection coefficient. 
Note that the measured response for larger plate size is off considerably 
since the far field condition is not met and the anticipated smaller return 
is obtained. These measured reflection coefficients were obtained from cal- 
ibrated scattered fields using a 5” square conducting plate as a reference 
standard. The error shown in Figure 2.4 can be minimized if it is cali- 
brated with a reference plate of the same size. All measurements employed 
vector background subtraction which entails the subtraction between mea- 
surements when a plate is present and absent. 
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Figure 2.2: Far field scattering for a plate with plane wave incidence. 
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Figure 2.3: Measured (solid) and calculated (dashed) reflection coefficients 
for a 5” square planar sample of lossy magnetic material. 
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Figure 2.4: Measured (solid) and calculated (dashed) reflection coefficients 
for a 12” square planar sample of lossy magnetic material. 
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Chapter 3 

Material Parameter Extraction 
Using Waveguides 


A common approach for parameter determination is to obtain reflection and 
transmission measurements for a fixed length of material in a waveguide 
fixture (coaxial or rectangular) with the fixture cross section completely 
filled. The parameters extracted from such an approach are subject to error 
from sample condition ( homogeneity, dimensions and shape distortions) 
and measurement hardware calibration. Examples of error magnitude will 
be presented for some of the errors. The error occured from an incompletely 
filled fixture cross section is presented in the following chapter. The basic 
equations used for such parameter extraction are presented here. 

The measurements are facilitated with automated broadband network 
analyzers providing Sn (reflection) and S 21 (transmission) characteristics 
for a sample length of d [5]. These measurements can be related to another 
set of reflection (r) and transmission (T) characteristics for a material 
junction semi-infinite length. The relationships are 

_ (1 -r 2 (e,/i))r(e,A*) (3.1) 

bn " (l-T’(£^)r 2 (e,/x)) 

and 
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(3.2) 


(1 - I^e./Ome,/.) 

3,1 (i-r>(e.M)r>(£,/i))' 

The previous expressions can be rewritten to express T and T in termB 
of Sn and Sn. These expressions are 


Sn -i- Sai — r 
i — (s n + s 21 )r 


(3.3) 


and 


r = K ± y/K~ 3 - 1 


(3.4) 


where 


_ s; 1 -s » 1 + i 

2S„ 


(3.5) 


The sign in Equation (3.4) is chosen to insure |T| < 1. 

These new characteristics, T and T, can be related to the material pa- 
rameters. The relationships for a coaxial waveguide fixture are 


T = 


(3.6) 


and 

(3. 

Note that t — t T t Q and \i = ji r /x 0 where the subscript notation refers to 
relative and free space values. 

The corresponding relationships for a rectangular waveguide fixture are 
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-i^Xi^tTJTrd 


( 3 . 8 ) 


T = e 


and 


r = 



( 3 . 9 ) 


where X. = yfl -(&) ! and X, = with *• bein « the tree 

space wavelength, c the speed of light and a is the width of the rectangular 

guide. 

The final expressions for the constitutive parameters using a coaxial 
waveguide are given by 



( 3 . 10 ) 


and 


£ r 



where 



and 



( 3 . 11 ) 


( 3 . 12 ) 


( 3 . 13 ) 


The final expressions for the constitutive parameters using a rectangular 
waveguide are given by 
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(3.14) 


_ i + r \ 0 
^ = i - r IF 0 


and 




f- + — 

U 1 (2a) 3 / 


(3.15) 


where 


1 

A 2 


2?r<i T 


(3.16) 


Note that Inj is modulo j2 , nn where n is the integer of The proper 
value of n can be estimated with the integer value of 


n = Int 


7<E' ~ 

27T 


(3.17) 


where <t>' is the derivative (slope) of the phase for the transmission coefficient 
T with respect to frequency and / is frequency. 

An important aspect in parameter determination is the error occurred 
from calibration and sample dimensions. The following examples are for 
coaxial and rectangular waveguide fixtures with sample parameters of fx r = 
(1.42, —j. 8) and e r = (7.8, — J.17) and a length of .05”. Figures 3.1 and 3.2 
illustrate error from the phase reference plane position for phase offsets in 
Sn corresponding to distances of .001” and .003”, respectively. The error 
is defined as 


% 


Error = 


Extracted Value - Exact Value 
Exact Value 


(3.18) 


The phase offset error (position) resulted in the same amount of error for 
both fixture types with the imaginary component of e being very sensitive 
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to the position of the reference plane. Similar results also occur for the n 
parameter. 

Figures 3.3 and 3.4 illustrate the error when the sample length is as- 
sumed to be .048” and .052”, respectively. The fixtures performed differ- 
ently in this case with the coaxial fixture yielding a lower error than the 
rectangular waveguide fixture. 

Other sources of error arise from other discontinuities occurring in the 
measurement system. The best solution so far has been to time gate the 
desired returns from the others which is readily performed with swept fre- 
quency measurements. It is also advisable to use a fixture with sufficient 
length to isolate the returns due to the fixture junctions from the sample 
return in the interior region of the fixture. 
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R(E) % ERROR 1(E) % ERROR 




Figure 3.1: Normalized error for e when 5u has a phase position er- 
ror of .001” (solid) and .003” (dashed). The sample parameters are 
fi r = (1.42, -j.8) and e r = (7.8,-j.l7) with d = .05”. Coaxial airline 
fixture. 
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Figure 3.3: Normalized error for e when Sn has a sample length er- 

ror of —.002” (solid) and +.002” (dashed). The sample parameters are 
/i r = (1.42, — j.8) and e r = (7.8, —j. 17) with d = .05”. Coaxial airline 
fixture. 
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R(E) % ERROR 1(E) % ERROR 




FREQUENCY IN GHZ 


Figure 3.4: Normalized error for e when 5n has a sample length er- 

ror of -.002” (solid) and +.002” (dashed). The sample parameters are 
fi r = (1.42, -;.8) and e r = (7.8,-j.l7) with d = .05”. Rectangular waveg- 
uide fixture. 
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Chapter 4 


Partially Filled Waveguides 


An accurate parameter extraction using a waveguide fixture with a cross 
section completely filled with material requires that no air gaps exist be- 
tween the material and the walls of the waveguide. It is believed that the 
tolerance of the dimensions of a material sample can be relaxed if the field 
solution has a defined air gap. Any slight changes due to sample dimension 
error or temperature expansion will not dramatically impact the measured 
quantities as significantly as if no air gap is to be present. 

The geometries under consideration are shown in Figure 4.1. Geometry 
B for the rectangular waveguide is more desirable than Geometry A since 
a small air gap along the short wall length is not as critical as along the 
longer wall length due to the electric field direction of the excitation field. 

The disadvantage for such a sample geometry presents itself in actu- 
ally acquiring the material parameters. The Su and Sji parameters have a 
complicated analytical representation which involve a moderate amount of 
numerical calculation to obtain. Hence, explicit expressions are not avail- 
able to obtain t and fi. The parameters e and fi have to be obtained through 
a numerical Newton-Raphson search. The search involves iterating on e and 
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fi until the following two functionals are driven to zero: 

771 rjcal cmeai 

*1 “ ^11 ” ^11 

jp Qcal emcaM 

*2 — ^21 “ °2l * 

The calculation of S n and S 2 i involves the modal expansion of the fields 
in all three regions and performing a modal matching at the material in- 
terfaces (z = 0 and z — d) for the geometry shown in Figure 4.2 [6]. The 
modal matching is essentially a three dimensional method of moments solu- 
tion using entire basis functions. The modal matching at the two interfaces 
were reduced to just one interface using a self-consistent approach. The 
required Sn and S 2i can be related to another set of scattering parame- 
ters for a semi-infinite geometry as shown in Figure 4.2. The relationships 
between these parameters are 

s„ = {[»„] + (.„][PMP) [[/] - MlflMliT 1 [*«]} (43) 

and 

Sn = [*„][P] [|J| - MPMPir 1 M- (4.4) 

The variables ay are the modal scattering matrices for the semi-infinite 
geometry and [P] is the modal propagation matrix corresponding to length 

d. 

I Field Expansions 

The required elements for the scattering matrices [sy] are obtained through 
a modal matching technique for the field expansions in both regions. The 


(4.1) 

(4.2) 
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tangential field matching yields expansion coefficients for the expansion 
modes. The expansion modes for the Geometry A in a rectangular fixture 
are given below. 

Region I 

E y = (e~ ri ‘ + aie r ‘*)$i + f) a n $ n e r "* 

n=2 


( 4 . 5 ) 


and 


H x = -Y x (e- r " - a 1 e r '*)* 1 + £ o n y n $ n e r -** 

n=2 

where 



Region II 


E y = ^b n ^ n e-^ 

n=l 


and 


H m = -Y d bnY in ^ n e-^ t 

n= 1 


( 4 . 6 ) 


( 4 . 7 ) 

( 4 . 8 ) 

( 4 . 9 ) 


( 4 . 10 ) 


( 4 . 11 ) 


where 


= C sin k m , n x for 0 < * < d, 

sinfe «»n( a - g) for d < x < a, 
7n = ^*i n bi\e T ^ T = fc* 0n — fc*, 

for 0 < * < d, 

Y 0 n = 7z% for d < x < a 


( 4 . 12 ) 
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and C is a normalization constant. 


Comparable field expansions are required for the other geometries but 
are not presented here for simplicity. The unknown coefficients are then 
obtained by solving for them with following expression 

[ Ex H dS = / Ex H dS (413) 

Jt= 0 - Jz=0+ 

through generating a system of equations from a method of moments ap- 
proach using entire bases testing functions. 

II Calculations and Measurements 

Software only for the rectangular Geometry A has been completed as of 
now. Figures 4.3 and 4.4 illustrate Sn and S 21 for t = .2 , .5 , .7 , and 
.9” in a X-band fixture (.4”x.9”) when fi r = (1.42, -j. 8), e r = (7.8, -j.17) 
and l = .3”. The variation of Sn and S 21 decreases rapidily as the guide 
becomes completely filled which indicates a small air gap to the short wall 
may not be that critical. 

Figures 4.5 and 4.6 illustrate Sn and S 21 for £r =2,5,8,11, and 14 in 
a X-band fixture (.4”x.9”) when fi r = (1.,-jO.), t = .1” and / = .5”. The 
variation of S n is quite sensitive for any loading and S 21 would show more 
variation if loss was present in the sample. 

Figure 4.7 illustrates the error in extracting the parameters assuming 
the cross section of the waveguide is completely filled but actually having 
an air gap. The air gap was .002” with fi T = (1.42, —j.8), e r = (7.8,— j. 17) 
and / = .3”. As anticipated, the error iB very small since the electric field 
in the air gap is parallel to a conducting boundary and hence is very small. 
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The error shown is comparable to the results when the fixture is completely 
filled giving an indication of the accuracy of the calculation. 

Figures 4.8 and 4.9 illustrate the comparison of the measured and calcu- 
lated Sn and S 21 parameters for a X-band waveguide shown in Figure 4.10. 
The sample dimensions were t = .082” and / = .9”. The calculated response 
used estimated parameters of \i T = (1.42,— j. 8) and e r = (7.8,— j. 17). The 
comparison is good and any differences can be attributed to the estimated 
parameters used in the calculations. The parameter estimates were “aver- 
age” values obtained from coaxial fixture measurements. 

Comparable results for the rectangular waveguide Geometry B have 
not been obtained yet. Figures 4.11 and 4.12 are measured results for 
the rectangular waveguide Geometry B for comparison purposes to the 
Geometry A results. 
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Figure 4.2: Waveguide regions. 
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Figure 4.3: Calculated S n for t = .2”, .5”,. 7”, and .9” in a X-band fixture 
(.4”x.9”) when [i T — (1.42, —j. 8), e T = (7.8,— j.17) and l = .3”. 
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Figure 4.4: Calculated Sji for t = .2”, .5”, .7”, and .9” in a X-band fixture 
(.4”x.9”) when fi r = (1.42, -j. 8), e r = (7.8,-j.l7) and / = .3”. 
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Figure 4.5: Calculated Sn for £ r = 2,5,8,11, and 14 in a X-band fixture 
(.4”x.9”) when fi r = (1., — jO.), / = .1” and / = .5”. 
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Figure 4.6: Calculated S 2 i for e r = 2,5,8,11, and 14 in a X-band fixture 
(.4”x.9”) when /i r = (1.,— jO.), t = .1” and / = .5”. 
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R(E) % ERROR 1(E) % ERROR 



Figure 4.7: Calculated error for a .002” air gap along the short wall for a 
rectangular guide with fi T = (1-42, — j.8), t T — (7.8, — j.17) and / = .3”. 
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FREQUENCY IN GHZ 


Figure 4.8: Measured and calculated Sn for a partially filled rectangular 
X-band guide with t = .082”. Geometry A. 
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Figure 4.9: Measured and calculated S 21 for a partially filled rectangular 
X-band guide with t = .082”. Geometry A. 
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Figure 4.10: Illustration of X-band fixture used for material measurements. 


ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPH 


34 



MAGNITUDE IN DB PHASE IN DEGREES 



Figure 4.11: Measured Sn for a partially filled rectangular X-band guide 
with t = .082”. Geometry B. 
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MAGNITUDE IN DB PHASE IN DEGREES 



FREQUENCY IN GHZ 


Figure 4.12: Measured S 21 for a partially filled rectangular X-band guide 
with t = .082”. Geometry B. 
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Chapter 5 
Conclusions 


The major focus for this work has been to examine alternate sample ge- 
ometries to relax the required tolerances on sample dimensions for elevated 
temperature measurements. The numerical effort for the partially filled 
fixture is considerable more than the completely filled fixture. The numer- 
ical solution for the other geometries are presently being generated. It is 
expected that an air gap for the other geometries will be more critical due 
to the gap orientation and electric field direction. The most appropriate 
partially filled geometry for parameter extraction will be found based upon 
numerical effort and final fixture considerations for elevated temperature 
measurements. 
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